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The strain and electric fields present in free-standing AlGaN/GaN slabs are examined theoretically within the 
framework of fully-coupled continuum elastic and dielectric models. Simultaneous solutions for the electric 
field and strain components are obtained by minimizing the electric enthalpy. We apply constraints appropriate 
to pseudomorphic semiconductor epitaxial layers and obtain closed-form analytic expressions that take into ac- 
count the wurtzite crystal anisotropy. It is shown that in the absence of free charges, the calculated strain and 
electric fields are substantially differently from those obtained using the standard model without electromechan- 
ical coupling. It is also shown, however, that when a two-dimensional electron gas is present at the AlGaN/GaN 
interface, a condition that is the basis for heterojunction field-effect transistors, the electromechanical coupling 
is screened and the decoupled model is once again a good approximation. Specific cases of these calculations 
corresponding to transistor and superlattice structures are discussed. 

PACS numbers: 73.21.Ac, 71.20.Nr, 85.30.Tv, 73.20.At, 73.61.Ey, 77.65.Ly 



I. INTRODUCTION 

In previous analyses of the strain in epitaxial layers of 
AlGaN/GaN, the electrical and mechanical properties of the 
crystal have been treated as though they are decoupled. Linear 
elastic theory is assumed to hold, and Hooke's law is invoked 
to obtain the relation between the in-plane and axial compo- 
nents of the strain tensor The standard model decouples the 
strain tensor from the electric field and enables the separation 
of the electric field and electronic eigenstate calculations from 
calculations of the strain field. Historically, there is a solid 
tradition of using this separability when studying heterostruc- 
tures of GaAs and associated alloys where the approximation 
that the strain field is negligibly affected by the electric field 
is validii 

Generally, it is known from thermodynamics^ that the elec- 
trical and mechanical properties of piezoelectric materials are 
coupled and a simultaneous treatment is called for when the 
electromechanical coupling is significant. This is the case 
when the piezoelectric response is large, as it is in AlGaN, es- 
pecially for high Al fractions. For example, large corrections 
of the electrostatic and elastic properties have been predicted 
for AlGaN/GaN transistor structures^ and for free-standing 
plates of AIN"^ when a fully-coupled model is used instead 
of a standard (uncoupled) one. 

In this paper, we calculate the strain field in a free-standing 
AlGaN/GaN slab using a fully-coupled model. A free- 
standing bi-layer slab is chosen as the model structure on 
which to develop the theory as it can serve as a building block 
for other, more complicated, structures. For instance, it can be 
used as a building block for superlattices and quantum wells. 
In addition, a heterojunction field-effect transistor (HFET) is 
just a special case of a bi-layer slab in the limit that the GaN 
layer is much thicker than the AlGaN layer It will be seen 



that the fully-coupled strain and electric fields for an HFET 
can be obtained in this asymptotic limit. 



Although the present model is still within the framework of 
continuum mechanics, it goes beyond the standard continuum 
elastic theory used in typical strain calculations in semicon- 
ductor materials. We obtain the total electric enthalpy for the 
bi-layer slab and minimize it subject to certain constraints to 
find both the strain and electric fields. Two constraints are 
used. One is that the two layers must share a common c-plane 
lattice constant after strain. This lattice constant is unknown at 
the outset and is deduced only after minimization. The other 
constraint is set up by the electrostatic potential, and hence the 
electric field, being forced to satisfy the Poisson equation sub- 
ject to the boundary conditions at the surface and the common 
interface. Both the spontaneous and piezoelectric polariza- 
tions are included. The effect on the strain and electric fields 
of a two-dimensional electron gas (2DEG) at the AlGaN/GaN 
interface, a situation that is essential for channel conduction 
in the HFET, is also investigated. The model produces ex- 
pressions for the strain tensor and the electric field in the two 
layers in closed form. 



This paper is organized as follows: In Sec. HI] the electric 
enthalpy for the AlGaN/GaN slab is minimized subject to the 
constraints derived from the pseudomorphic strain condition 
and from the Poisson equation. The strain tensor and the elec- 
tric field are derived in closed form. In Sec. |lll| the calcu- 
lated results from the fully-coupled model are contrasted with 
those from the uncoupled model. The effect on the strain and 
electric fields of screening from the 2DEG is discussed. The 
results are summarized in Sec. lIVI 
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II. MODEL DESCRIPTION 

In the standard strain theory for semiconductor materials, 
the formal way to calculate the strain field in a generalized 
strain problem is to minimize the Helmholtz free energy^ 

F ^ ^Cijknijjki, (1) 

in which Cijki is the fourth-ranked elastic stiffness tensor, 7^ 
is the strain tensor, and the indices i, j, k, and / run over 
the Cartesian coordinates x, y, and z. Summation over re- 
peated indices is implied throughout. The minimization will, 
of course, be subject to certain constraints brought about by 
the boundary conditions at the surfaces and between adjacent 
materials, as well as any external forces. 

In a wide range of semiconductor problems, however, one 
is concerned with calculating the strain in epitaxial layers 
grown, at least in principle, pseudomorphically on a thick 
buffer layer Typically, the thick layer belongs to, or at least 
closely resembles, the same crystallographic point group as 
the epilayers. In this case, a minimization of Eq. Q is unnec- 
essary, since we can take full advantage of the boundary con- 
dition for a free surface which requires that the components of 
the stress tensor parallel to an outward normal to the surface 
be zero. For a surface normal to the z-axis, this means that 
o'iz = 0, where is the stress tensor In addition, since the 
lateral extent of the layers is far greater than their thickness, 
a one-dimensional (1-D) approximation can be used, wherein 
•jij ~0 for i^ j. This condition is equivalent to ignoring the ef- 
fects of bowing on the local lattice displacement. Under these 
conditions, the result 

CTzz = = Cxxzzhxx + lyy) + Czzzz^zz, (2) 

is obtained for wurtzite epilayers oriented in the [0001] di- 
rection. Here the only unknown is 7^^, because a stan- 
dard approach is to assume that the buffer layer is unstrained 
and that the epilayer assumes the c-plane lattice constant of 
the buffer, effectively fixing ^^x and ^yy. Clearly, when 
^xx — Tyy 

= {a~ao)/ao, Eq. 12 yields the usual resuh for the 
Poisson effect in the case of biaxial strain in an epilayer which 

is ^zzl^xx '^Cxxzz I C zzzz- 

The situation for the free-standing bi-layer slab depicted in 
Fig- [Din which the two layers are of comparable thicknesses 
is more complicated, since none of the strain components is 
known beforehand. To obtain the strain, at least within the 
standard model, the total mechanical strain energy for the slab 
must be minimized with respect to 7^ subject to the con- 
straints along the interface. There is a further complication 
in piezoelectric materials such as AlGaN and GaN wherein 
the electrical and mechanical properties are coupled through 
the piezoelectric coefficient tensor Equation Q is no longer a 
suitable energy functional for calculating the strain field, since 
it does not include the electromechanical coupling. Instead, 
we begin the derivation with the electric enthalpy H given 
by^ 

iJ = L/-E-D, (3) 



where E and D are the electric field and electric displacement, 
respectively, and U is the total internal energy given by 

^ = 'i^Cijknijlki + ]^eijE.iEj, (4) 

in which is the tensor form of the electric permittivity. The 
electric displacement field components are given by 

Di = e.jkljk + £-ijEj + PI, (5) 

in which e^fe is the piezoelectric coefficient tensor and Pf 
is the spontaneous polarization.** The first term in Eq. (jSj is 
recognized as the piezoelectric polarization. The spontaneous 
polarization is in the direction N ^ Ga along the c-axis bond. 
After substitution into Eq. (|3|l, the enthalpy in its final form is 

P[ = -^Cijknij^ki — BijkEi^jk — -EijEiEj — EiP-. (6) 

This expression includes the full electromechanical coupling 
as well as the spontaneous polarization. Differentiating the 
enthalpy with respect to the strain tensor gives the fully- 
coupled equation of state 

fij = CijkHki — SkijEk, (7) 

for piezoelectric materials. To obtain the strain and electric 
fields within a fully-coupled model, Eq. (|6} must be mini- 
mized with respect to both 7^ and Ei subject to the con- 
straints to be discussed shortly for the bi-layer slab of Fig. 
□ 

Using the Voigt notation^ and expanding Eq. (|6|l, the en- 
thalpy for a wurtzite crystal with the [0001] axis as the princi- 
pal axis can be written as 

H = ^Cii(7^^ + -fly) + ^Csslzz + Cl2jxxlvy + 
Cisjzzilxx + lyy) + IC^^ijl^ + 7yJ + 
(Cii - C-Y2)liy - EzV'ilinxx + lyy) + 6337^^] " 

ex^{Exixz + Ey-fyz) - EzP' - \^e{El + E^ + El), (8) 

where e is the electric permittivity assumed here, for simplic- 
ity, to be a scalar The enthalpy for the slab is given by 

ta + U 

where ta and ti, are the thicknesses of the two layers and "a" 
denotes the AlGaN layer and "b" the GaN layer It is noted 
that the elastic and piezoelectric coefficients and the electric 
permittivity can be substantially different in the two layers. 
These variations are taken into account in Eq. (|5Jl. Because of 
the rotational symmetry of the slab, jyy = jxx and •jyz — ixz 
in each layer The shear terms in the strain tensor, 7jj for 
i ^ j, turn out to be zero because of the assumption that the 
strain is piecewise homogeneous. A more realistic assumption 
would be to expect the strain to diminish away from the inter- 
face, resulting in a non-zero •jxz that would vary with position 
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along the c axis. This situation would be manifested in a bow- 
ing of the slab. Such inhomogeneous strains, however, would 
require a three-dimensional (3-D) numerical calculation that 
is beyond the scope of the present work. 

One of the minimization constraints on Eq. (|9j is that the 
two layers share the same c-plane lattice constant which is 
unknown at the outset. This condition is expressed as 



and 



(1 +7Daa = (1 +7Dab: 



(10) 



where aa and ab are the unstrained c-plane lattice of the Al- 
GaN and GaN layers, respectively. Equation dlOt presumes 
ideal growth conditions. Partial relaxation due to dislocations 
will be the subject of future investigations. 

The other constraint is the relationship between the electric 
and strain fields. This relationship is established in the present 
work by solving the 1-D Poisson equation given by 

d d(f) d .^^ ^ , ^, , 

^e— = ttI-P + 26317x2: + 6337^^) + 6071200(2; - ta), 
oz az oz 

(11) 

where (f) is the electrostatic potential and n2D is the sheet con- 
centration of an ideal 2DEG modeled as a (5-function at the 
AlGaN/GaN interface. Equation ( II H is solved analytically 
for (f) subject to the boundary conditions = at z = and 
z ^ (ta + tb) and also to the continuity of (j) and the electric 
displacement across the interface. The general solution of Eq. 
il 1> in each layer is given by 

P' 263172;^ + 6337^^ A 

(j) = — zH zH z + B, (12) 

e e e 

where A and B are unknown constants. It is evident that there 
are four unknowns, two in each layer The S's are eliminated 
by enforcing the boundary conditions </) = at z = and 
z = ta + tb- This condition at the two surfaces presumes that 
the polarization charge of the materials is terminated by exter- 
nal charges and that the applied bias voltage is zero. Further, it 
ensures that the electric field in the free space above and below 
the slab is zero, a reasonable and physically plausible result. 
The relationship between the two A's is established from the 
continuity of the electric displacement and is obtained by in- 
tegrating Eq. (II It across the interface (see Fig.^ as follows: 



dz 



PT"_ + (26317... + 6337..)! - + 6o7l2D- (13) 



From Eqs. M2\ and (I13> . the relation between the A's in the 

respective layers is given by 



AlGaN 



60"2D- 



(14) 



From Eq. ( I14> and from the continuity of at z = ta, we can 
solve for all of the remaining unknowns. 

With the electrostatic potential now known in terms of the, 
as yet, undetermined strain tensor, the electric fields in the two 
layers are then given by 



Et, = tb 



- 2e^i7^, + 2e^i7; 



'xx ^Salzz + 6337zz 



ta£b + tb^a 



(15) 



P"-' - 26^i7L + 2e|i7L - el,jl 

ta^b + tb^a 



6337ZZ 



El = -ta 

tnth + thtn 

(16) 

where P"'=' = P***' - P'**''' + eonio- Equation ^ can be 
readily minimized with respect to 7^ and subject to the 
constraints expressed in Eqs. dlOL M5\ . and M6\ using com- 
mercial symbolic-mathematical programs. For convenience, 
we have also formulated the problem in the form of a matrix 
equation of the form 



All A12 Ai3 
A21 A22 A23 
,A3i A32 A 



33 




(17) 



in which the matrix elements are obtained by differentiating 
the energy functional with respect to the unknown strain com- 
ponents. The constraints are already built into Eq. ( I17> by 
eliminating 7^^., El, and E^ from Eq. (|9}. Equation (I17> can 
be solved either symbolically or numerically. We have pro- 
vided the symbolic solution as an auxiliary itemi^ Appendix 
|X]lists the matrix elements. 

The fully-coupled model described above reproduces the 
strain tensor from the standard model in the limit that the 
piezoelectric stress tensor is set to zero. The results are as 
follows: 



^xx I e 



ijk- 



and 



for the AlGaN layer and 



and 



7' I 

/ XX 



7" I 



tbO-aio-b - 


- aa)C^3 




taOlC^^R 


^ tbUaiab 


— aa)C'i3 


tbalC^s - 


~ taCl^C^^R 


taO,b{o.b - 


- aa)C'33P 


tbalC^,, 4 


- taalC\^R 


^ taab{ab - 


0-a)C\^R 




taalCl^R 



(18) 



(19) 



(20) 



(21) 



for the GaN layer where the constant factor R is defined as 



R 



(Cii + Ci2)C33 — 2C; 



-"13 



(Cll + Cl2)C3; 



33 



2C?3 



(22) 



It is further seen from Eqs. ( I18> -( I2U that in the limit tb ^ 
ta, the well-known uncoupled result for the axial strain in an 
epitaxial AlGaN layer pseudomorphically strained on a thick 
GaN buffer is recovered: 



7; 



= -2 



!xx\tb—^oo 

'-'33 



(23) 



In this limit, the strain in the GaN layer is reduced to zero, at 
least ideally, as the thick GaN layer serves as a fixed anchor 
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The uncoupled electric field in the AlGaN layer in this limit 
then becomes 



El 



pnet 



'^33 "-^13 



^xx \t 



(24) 



We can also extract the fully-coupled results for the strain 
and electric fields in the limit tf, ^ ta, i.e. the usual HFET 
configuration, and compare them with the uncoupled results. 
In this case, fully-coupled axial strain in the AlGaN layer is 
given by 



^zz too 

2e33(e|3C23 — e'^iC^s) 



C: 



Ix 



33 



^33(60(733 + e' 



-33 / 



^xx 



— too 



+ 



a pnet 



33 



-33 



and the fully-coupled electric field by 



El 



C^3^'"' 

f r"^ -I- pa 2 
taL'33 ^ ^33 



a /^a _ a /^a 
^31^-^33 ^33'-^13 

p r'a 4- pa 2 
ta^-^33 ^33 



^xx\th—>c 



(26) 



In Eqs. ( I23t . ( I24> . ( I25> and ( I26> . 7^^ takes the asymptotic 
limit (cb — aa)/aa as — > 00. The first term in Eq. i25\ is 
the result from the standard model. The other two terms are 
due to electromechanical coupling. 

It is worth pointing out that the foregoing results also apply 
to a free-standing AlGaN/GaN superlattice in which the slab 
of Fig.^forms a period of the superlattice. This outcome fol- 
lows from enforcing the pseudomorphic boundary condition 
throughout the entire cross-section of the slab, i.e. for all z. In 
turn, this result is a consequence of the simplifying assump- 
tion made at the outset that 7^ = for i ^ j. In addition, 
the electrostatic potential will be subject to periodic boundary 
conditions at z = and z = ta + tf,. We can assume without 
loss of generality that the potential in the superlattice case can 
be set to zero at the two ends, as we have done in the present 
case. A non-zero potential at the two ends in the periodic sys- 
tem simply represents a rigid shift of the potential distribution 
function and does not change the electric field. Thus a super- 
lattice formed from a stack of slabs discussed in the present 
work will have identical strains and electric fields within each 
period as obtained for our model structure. 



III. RESULTS AND DISCUSSION 

Table |I] shows the material parameters ^^•'■Siiiii^ii^iiiii^ used 
in the calculations. The signs of the polarization parameters 
are defined in relation to the [0001] direction: a negative sign 
means that the vector is in the [0001] direction. 

For the model structure of Fig. ^ we choose ta = 300 A 
and tb = 500 A. The 300 A AlGaN layer is typical of most 
high-power HFETs. It is recognized that the GaN layer used 



in the model structure is much thinner than that permitted by 
current technology. For example, using a laser lift-off pro- 
cess and growth via metal-organic chemical vapor deposition 
(MOCVD), free-standing layers of nitrides can be produced 
successfully only for relatively thick layers'* in the vicinity of 
5 /im. The free-standing layers are even thicker for epilayers 
grown by hydride vapor phase epitaxy (HVPE), reaching be- 
tween 250 - 300 /im.'^ '* The fully-coupled model presented 
herein is quite general and, as shown in Sec. |ll] can readily 
reproduce the results for thick GaN. Figure |2 shows the cal- 
culated strain tensor as a function of the Al fraction using the 
fully-coupled model for a model structure that is depleted, i.e. 
«2D = 0. Following standard convention, a positive sign in- 
dicates dilation and a negative sign contraction. Unlike the 
situation of a AlGaN layer on a semi-infinite GaN buffer, both 
layers are now strained, with the in-plane strain tensile in the 
AlGaN layer and compressive in the GaN layer 

Next we examine the deviation between the fully-coupled 
and uncoupled models for a depleted slab. This deviation A 
is defined as 



a(b) a(b) uncoupled 

. a(b) "fa ~ 7m 

ii ~ a(b) 



(27) 



where i is x or z, the superscript "a" or "b" refers to a par- 
ticular layer, and 7^,; is given by the solution to Eq. illi for 
the coupled case and by Eqs. ( I18> - i2l\ for the uncoupled 
case. The calculated results are shown in Figs. |3] and |3 for 
the in-plane and out-of-plane strains, respectively. It is seen 
from Fig.l^that the deviation for 'y^x is quite small. A possi- 
ble reason is that the electromechanical coupling of '^^x into 
the equation of state [Eq. Q] occurs through the piezoelectric 
coefficient 631 which is much smaller than 633, as seen from 
TableU The electromechanical coupling of j^z into the equa- 
tion of state is through 633 which is quite large, especially for 
AIN. It is clear that the error in the calculated strain along the 
growth axis is quite significant, reaching 30 % in the AlGaN 
layer for a mole fraction of 0.3. This error will, in turn, have a 
significant impact on the calculation of electronic and optical 
properties that depend on accurate knowledge of the strain. 
Some examples of such quantities include the eigenstates and 
the complex dielectric function of the slab. 

So far it has been assumed that the slab is depleted with the 
only charge present being a fixed space charge P given by 



p ^ ps (b) _ ps (a) _|_ pp (b) _ pP (a) 



(28) 



where and PP are the piezoelectric polarizations 

PP(^> =2e^i7L + e^37L, (29) 



and 



pp*' =2e^3i7L 



^33! XX- 



(30) 



For a cation-faced structure, the standard growth orientation 
for most HFETs, P is positive. '"^ Depending on the growth and 
surface conditions, donors may be present and may contribute 
free electrons to the interface. Under certain conditions, 
the space charge may be partially screened by a 2DEG. For 
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instance, using our own self-consistent Schrodinger-Poisson 
model,-" we find that the 2DEG is almost 90 % of the mag- 
nitude of the fixed space charge. We assume, conservatively, 
that ri2D = 0.8P/eo and examine the effect of free-electron 
screening on the electromechanical coupling. Figures |5] and 
|6lshow and 7^^, respectively, calculated using the stan- 
dard and coupled models. The in-plane strains are not shown, 
since the effect of electromechanical coupling on the in-plane 
strain components is weak. It is evident from the calculated 
strains that the effect of coupling is to reduce the magnitude 
of the strains relative to those of the standard model, espe- 
cially when the structure is depleted. This result is not surpris- 
ing, since we would expect the piezoelectric fields to induce 
forces to oppose any applied forces such as those due to the 
pseudomorphic interface condition. Another result that is evi- 
dent from Figs.|5land|6lis that the 2DEG screening brings the 
strains of the coupled model closer to those of the standard 
model. As will be seen shortly, the effect of screening is to 
reduce the built-in electric fields. This, in turn, will reduce the 
electromechanical coupling, as seen from Eq. Q. 

Figures and |8] show the calculated electric fields in the 
AlGaN and GaN layers, respectively, using the standard and 
coupled models. These fields help to explain, in part, the large 
deviation of the standard strain fields from their coupled coun- 
terparts. As an example, for a Al fraction of 0.3, the electric 
field in the AlGaN layer is about 2 MV/cm, provided that the 
structure is depleted. Such a large field, if present, will give 
rise to a strong electromechanical coupling and, therefore, a 
large discrepancy between the standard and coupled models. 
When there is a 2DEG, however, the large space charge that 
gives rise to the electric field is partially neutralized, depend- 
ing on the magnitude of the 2DEG. For n2D = 0.8P, the field 
in the AlGaN layer is reduced to about 0.4 MV/cm for a Al 
fraction of 0.3. Evidence for the smaller field is seen in re- 
cent photoreflectance (PR) measurements2i22. on AlGaN/GaN 
heterostructures in which barrier fields in the vicinity of 0.3 
MV/cm were reported for 250 A thick barriers with Al frac- 
tions in the range 0.06 -0.1. An electric field of opposite sign 
occurs in the GaN layer with similar trends between the cou- 
pled and standard models. Its magnitude is decreased with the 
thickness of the GaN layer. It is recognized that the calcu- 
lation of the electric field in the GaN layer is more involved 
than the classical model presented here allows. The formation 
of the notch in the conduction band edge in which the 2DEG 
resides is reproducible only within the scope of a quantum- 
mechanical calculation. In spite of approximations used, the 
calculated fields shown here appear to be quite plausible. 



components and electric field in each layer are obtained. The 
results for the bi-layer system are also appropriate for free- 
standing binary superlattices with negligible bowing. Large 
discrepancies between the standard and coupled models are 
found for depleted structures due to large built-in electric 
fields. The electric fields are reduced when the polarization- 
induced space charge is screened by a 2DEG. As a conse- 
quence, the discrepancy between the standard and coupled 
models is reduced when a 2DEG is present. 
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APPENDIX A: MATRIX SOLUTION OF THE STRAIN 
FIELD 

The matrix elements required to solve Eq. illl are given in 
this appendix. As before, "a" refers to the AlGaN layer and 
"b" to the GaN layer. 

All = taO-liCli + Cl2){ta<^a + h^b) 

+ hal{taeb + Uea)(C\^+C\2) 

+ 2tah{el^aa - e^la^,)^ (Al) 

A12 = taalCl^{ta€b + hf-a) 

+ tatbabelr^{e\^ab - e\^aa), (A2) 

^13 = tbaaa,bC\^{tatb + h^a) 

- tatbabe\^{e%^ai, - e\,^aa), (A3) 

A21 = 2taabCl^{ta€b + tb^a) 

+ 2tatbel3{eliab-eliaa), (A4) 



A22 = taabC^sita^b + tb(-a) + tahabel^'^ , (A5) 



IV. SUMMARY AND CONCLUSIONS 

In summary, a fully-coupled electromechanical model has 
been presented for the strain and electric fields in a free- 
standing bi-layer AlGaN/GaN slab. The electric enthalpy of 
the slab, with contributions from the piezoelectric and sponta- 
neous polarizations, is minimized subject to the constraints 
imposed by the pseudomorphic interface condition and the 
Poisson equation. Closed-form expressions for the strain 



A23 = -tahabel^el^, (A6) 

^31 = 2tatbe^^{e\iaa - e^iab) 

+ 2tbaaC\3{taeb + tbea), (A7) 

A32 = -taibOfceggegg, (A8) 
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^33 — 4af)C'33(iae& + tb^a) + taha-be^s 



(A9) 



Bi = tatb{e%-^ab ~ egifla) 

+ aatb{ab - aa){C\^ + C\2){taeb + hea), (AlO) 



B2 = 2^0^663^633(00 — flfc) 

+ tahabel^iP'^'''' - F^(^) + 60^20), (All) 

and 

B3 = -tatbabel^iP''^'^ - P'^'^ + 6071215) 

+ 2{ab — ao)[iaibe3i633 + tbC^sitaeb + hea)]. (A12) 
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TABLE I: Strain-related material parameters used in the present 
model. The elastic stiffness constants are in units of GPa and the 
piezoelectric coefficients and spontaneous polarization in units of 
C/m^ 

Material Cn C12 C13 C33 a (A) 631 633 P" e/eo 
396° 137° 108° 373" 3.112* -0.58" 1.5; 
367° 135° 103° 405° 3.189' -0.36'' 1" 
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FIG. 1 : Cross-section of a model cation-faced AlGaN/GaN bi-layer 
slab showing the direction of the piezoelectric and spontaneous 
polarization vectors in relation to the z-axis. Both layers in the 
slab are under strain as shown, ta and tb are the thicknesses of the 
AlGaN and GaN layers, respectively. 
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FIG. 2: Calculated strain for the structure of Fig.Qas a function of 
the Al fraction with ta = 300 A, tt = 500 A, and 7120 ~ 0, using the 
fully-coupled model. 
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FIG. 3: Deviation of the uncoupled from the fully-coupled in-plane 
strain in the bi-layer slab with mo = 0. 
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FIG. 4: Deviation of the uncoupled from the fully-coupled out-of- 
plane strain in the bi-layer slab with n2D = 0. 
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FIG. 5: Calculated out-of-plane strain in the AlGaN layer for the 
standard and coupled models. Two coupled cases are shown, one 
without free electrons and one with free electrons. 
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FIG. 6: Calculated out-of-plane strain in the GaN layer for the stan- 
dard and coupled models. Two coupled cases are shown, one without 
free electrons and one with free electrons. 
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FIG. 7: Calculated electric field in the AlGaN layer for the standard 
and fully-coupled models showing the effect of free-electron screen- 
ing in each case. 
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FIG. 8: Calculated electric field in the GaN layer for the standard and 
fully-coupled models showing the effect of free-electron screening in 
each case. 



